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We comment on the problem of energy translation invariance of probabil- 
ity distribution and present some observations. It is shown that a probability 
distribution can be invariant in the thermodynamic limit if there is no long 
term interaction or correlation and no relativistic effect. So this invariance 
should not be considered as a universal theoretical property. Some peculiari- 
ties within the invariant g-exponential distribution reveal that the connection 
of the current nonextensive statistical mechanics to thermodynamics might 
be disturbed by this invariance. 
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I. INTRODUCTION 

The nonextensive statistical mechanics (NSM) [jT]-|7j is a generalization of 
the Boltzmann-Gibbs statistics (BGS) intended to study complex nonlinear 
systems having long-range correlations and fractal or chaotic state space 
for which BGS is no more valid. Due to the generalized character of this 
theory, the scientists of this field have been brought back to reflections about 
many fundamental aspects of physics theory. Can we keep the additivity 
of, e.g., entropy, energy and volume for complex systems? What are the 
possible relation between the probabilities of interacting systems? What is 
thermodynamic equilibrium in this case? How to define it for systems of same 
nonextensive nature and for systems of different nonextensive nature? Is the 
classical probability definition still valid in chaotic or fractal phase space? 
Does the completeness of information about systems containing independent 
parts still apply to systems with correlated or overlapped parts ... ? All these 
questions concerning the foundation of NSM are far from being answered and 
it would be too soon to conclude. A complete comprehension of these aspects 
is necessary to finally determine the validity limites of NSM and to know the 
kind of systems to which NSM should be apphed. 

One of the most worrying problems of the NSM distribution given by 
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P^ = ^[l-{l-q)^3e^]^^ (1) 

is that it is not invariant under uniform translation of the energy spectra 

Only when q = I, Eq.([I|) becomes pi = ^e~^^' and the invariance can 
be recovered. This g-exponential distribution can be given by maximizing 
Tsallis entropy^ with the unnormalized expectation []2| 

w 

U = Y.e^pl (2) 
i 

and the normalization Y.^ Pi = 1 from which the partition function Z is 
calculated. Eq.(|I|) can also be found, with a different partition function, 
by maximizing Tsallis entropy with Eq. (0) and an anomalous normalization 
J^f = 10- can be identified to the inverse temperature within a non- 
additive energy scenario [0-050] or with an additive energy approximation 

MB)- 

From the power law distribution Eq.([I|), it is obvious that, if we replace 
by e^ + C where C is constant, we find]3j(ei+C) ^ Pii^i)^ excepted that q = 1. 
If Pi varies with C, then the variance of all the thermodynamic functions will 
be disturbed by C. For example, for the energy U, we have 

U{e + C) = + CK(e. + C) (3) 

i 

^ E(e^ + C)[l-(l-gMe, + C)]T^ 



_ -1 dZ{e^C) 
~ Zi{e + C) dp 
_ d Z^-'jje + C) - 1 

where U = ~-^ ^\lq^ is the unnormalized expectation without energy trans- 
lation (C = 0). This theoretical feature seems difficult to accept according 
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to some scientists Later, a normalized expectation E of energy was 

proposed on the basis of the so called escort probability 



E = hili^. (4) 



which, introduced into the maximization of Tsallis entropy as a constraint 
replacing Eq.(0), allows to obtain 



-9 



p. = [1 - (1 - -/Me. - 
Z 



with 



w 



Z = Y.[^-{l-q)[5{e,-E)]-.. (6) 



The distribution Eq.(|5|) is invariant under energy translation because, accord- 
ing to reference 0], if we add a constant to e^, we have the same constant 
added to E. Recently, it is argued [jT^ that, for a correct and successful 



maximization of Tsallis entropy leading to Eq.(|^), the expectation Eq.(g) 
as a constraint should be replaced by the unnormalized expectation Eq.(0) 
minus the expectation given by Eq.(g), i.e. U — Y^iP^E = Y.iP^{ei — E). 

In this paper, we will firstly discuss different conditions for the validity of 
energy translation invariance of probability distributions. Then we present 
some observations about the relevant problems or theoretical peculiarities 
within invariant NSM. The probability invariance with a generalized energy 
shift consistent with NSM framework is discussed. 



II. ENERGY TRANSLATION INVARIANCE OF MAXWELL-BOLTZMANN 

DISTRIBUTION 



As well known, BGS distribution is invariant under energy translation. 
This property conforms with the intuition that the particle distributions in 
a container will not be disturbed when, for example, you raise the container 
to a higher point or move it at a different speed. This invariance of BGS 
is directly related to that of classical mechanics which leaves the interaction 
potential completely arbitrary and so is invariant if we add a constant to 
the total energy. In addition, the exponential form of BGS allows following 
invariant property : 
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1 e"^^ e~^^' (.-Pie^+C) 

pie.) = -e-^- = ^^-^ = ^-^^^ = p{e. + C) (7) 

So Pi will not be changed if we add a constant to the energy spectrum e^. 
The partition function will be changed as follows : 

Z{e + C) = j: e-^(^^+^) = e-^^ E = e-^^^(e) (8) 
j j 

which gives following internal energy U : 

U{e + C) = -Q^lriZie + C) = -—lnZ{e) + C = /7(e) + C. (9) 

It is instructive to see the role of thermodynamic limit in this invariance 
when we study the velocity distribution of ideal gas. Suppose that Vi is the 
velocity of the l^^ particle of mass m with respect to the container or to the 
center of gravity of the gas. The total kinetic energy of the gas is given by 

e = T.lrnvl (10) 
I ^ 

Now we are in another system of reference in which the container moves at 
a velocity vq, we have 

e' = E ^^(^0 + vi)^ = ^m(E + E ^0 + 2 E ^0 • ^/)- (H) 
I ^ ^11 I 

In the thermodynamic limit, the third term at the right hand side of Eq.(|TT]) 
is null because YlfLi vi = if N ^ oo. We then have e' = |m E/ vf^^mNif^ = 
e+^Mi)Q where M is the total mass of the gas. Put this in Eq.(|7D in replacing 
the summation by an integration in vi, we find that the distribution law does 
not change with vq. 

Briefly, according to the above discussion, we have at least three conditions 
on the invariance of the distribution. 

1. The mechanics basis must be non relativist, i.e. an arbitrary constant 
can be added to the energy. 

2. Thermodynamic limit holds, i.e. oo. 

3. The distribution law must be exponential. 
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It is well known that the third condition needs the second one and, in 
addition, the harsh assumptions that the interaction or correlation in the 
system is negligible or of short term so that the different parts of the system 
under consideration are independent and additive. In view of the above men- 
tioned conditions, it is difficult to say that the energy translation invariance 
of probability distribution is an universal theoretical property. Indeed, if we 
consider for example the relativistic effect, we can no more add constant to 
the energy e due to its relation with the total mass M given by 

e = Mc^ (12) 

where c is the light speed, e is not arbitrary because M can not be changed 
arbitrarily. If we suppose a system composed of many elements of mass mi 
and energy {i = 1, 2, 3, ...) in interaction with total potential energy M 
will be given by : 

M = Y.m, + V/c^ (13) 

i 

or 

Mc^ = Y.mc^ + v = Y.^^ + y (14) 

i i 

with Ci = rriiC^ for i^^ element. It is obvious that M, and V cannot be 
changed, or the variance of the theory would be perturbed |T3 . 



On the other hand, the relativistic kinetic energy of particles 
[mc^ / \Jl — vf /c^) is not linear function of vf. So what we did in Eq.(|TT]) 
leading to e' = e + mv\l2 does not hold for relativistic gas. This means that 
the relativistic velocity distribution will change with even with exponential 
distributions. 

Now if we consider complex systems with long range interactions which are 
not limited between the walls of the containers, the things are more complex 
than with short range interactions because the entropy and energy can be 
nonadditive so that the exponential distribution does not exist any more. So 
the invariant distribution is only a theoretical property in special cases. It 
should not be claimed for general cases. It is worth emphasizing that, up 
to now, all the successful apphcations of the distribution Eq.(|I]) are never 
disturbed by the fact that it is not invariant with energy translation. On the 
contrary, the invariant distributions like Eq.(|5|) were shown to present some 
serious theoretical difficulties, as discussed below. 
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III. ENERGY TRANSLATION INVARIANCE OF NONEXTENSIVE 

DISTRIBUTION 



Nevertheless, there is up to now no any exphcit reason to completely 
exclude invariant nonextensive distribution which is recently discussed in a 
general way and claimed to be an universal property verified by any thermo- 
statistics with linear or normalized expectation, whatever the entropy form 
We have shown in the previous section that this conclusion is not true. 
But a question remains open : is it possible to obtain the invariance with a 
distribution function like Eq.(^? In this section, we would like to show some 
observations about Eq.(^ and discuss some of its theoretical peculiarities. 

If we want a distribution pi to have energy translation invariance, we have 
to first of all define an invariant normalized expectation satisfying U{e + 
C) = [/(e) + C . This is possible with the expectation Eq.(|4D (or Eq.(@) plus 
the anomalous normalization E^J^f = 1 [j51-|7[]) under the assumption that 
Pi{e^ + C) = Pt{ei). That is 

E{e + C) = + C)p]{e, + C)/ Y.p'i{e^ + C) (15) 

i i 

= E + C 



Then the method of [T^ can be adapted to introduce J2f pjE as the in- 



variance constraint into the following functional 



w w 



A = -aY.Pr- PEple. ^iJ:p1E (16) 

^ Q i i i 



Let 1^ = 0, we obtain the following distribution : 

z 

For Pi to be invariant under energy translation, considering Eq. ([T5|) , we have 
to set /3 = 7, which leads to the invariant distribution Eq.l 



IV. WHAT IS WRONG WITH THE INVARIANT DISTRIBUTION? 

From the invariant distribution Eq.(5]), we easily show that 
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Z'-'EP^ = Erf[l - (1 - - E)] (18) 

i i 

Due to the invariant factor (e^ — E), Eq.([T^) leads to 

T.pI = z'-'^ (19) 

i 

and 

Z = U^-{l-q)P{e.-E)]T^.. (20) 

i 

or 

E[l - (1 - q)P{e^ - E)]Th = ^[1 - (1 - q)p{e, - E)]t^. (21) 

i i 

Eqs. (|raD and Eq. are obviously basic relations of the invariant theory 
and should hold for arbitrary value of (3 and e^. We will see that most of the 
problems encountered below are intrinsically related to these two equalities. 

1. Let us begin by addressing the problem of the calculation of nonextensive 
term or correlation in energy (or any other quantities of interest) with 
Eq.®. 

Suppose an isolated system C composed of two subsystems A and B in 
thermal equilibrium. It was shown [fH|J7| that, with Tsallis entropy, the 
equilibrium condition yields for even interacting subsystems the follow- 
ing product probability law 

p,,{C) = p,{A)p,{B) (22) 
and entropy pseudoadditivity 

S{A + B) = S{A) + S{B) + ^^S{A)S{B). (23) 

For NSM with g-exponential distribution, the product probability means 
that A and B are correlated and should give the pseudoadditivity as- 
sociated with the quantity of interest. But from Eqs. (|^) and (|22|), we 
straightforwardly obtain : 

e,j{A + B)- E{A + B) = [e,{A) - E{A)] + [e,(5) - E{B)] (24) 

+ {q-l)P[e,iA)-EiA)][e,iB)-E{B)]. 
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Without additional hypothesis, this equahty does not lead to any expHcit 
relation between the total energy and the energies of the subsystems A 
and 5, which is absolutely necessary for defining temperature, pressure 
and chemical potential as the measures of equilibrium. 

Recently, some authors [|TD|,|TT[ proposed neglecting the nonextensive 
term in energy and writing 

e,j{A + B) = e,{A) + ej{B) (25) 

and 

E{A + B) = E{A) + E{B). (26) 

We have seen that this additive energy approximation allows to reconcile 
the invariant theory with the zeroth law of thermodynamic [[TDIjlIIi 
establish a generalized thermodynamics. But it should not be forgotten 
that this is only an approximate approach. In other words, the invariant 
theory does not have vigorously defined temperature. This is one of its 
intrinsic flaws. 

2. Now we address a mathematical problem related to Eqs. ([T^ ) and (|2T|). 
From Eq. ([rap, Tsallis entropy can be recast as 

Z^"'^ - 1 

S = k- -. (27) 

1-q 



Then we calculate the following derivative 

dS k dZ 



dE Zi dE 
First we take the Z given by Eq. (g), we obtain 



(28) 



rl h rl w , 

U^-{l-qme^-E)]T-. (29) 



dE ZUE 

= gE[l-(l-9)/?(e.-^)]^ 

= kpz^-'^. 

But considering Eq. (|2T|) , we can also take the Z of Eq. (|2tl|), this time 
we obtain 
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dS qkB ,^ , ^, ,.xn29-i 

^ = ^ E[l - (1 - - ^)] — . (30) 

Comparing Eq.(|3D]) to Eq.(|29|), we get 

Z = qt[l - {I - q)(5{e, - E)f^ . (31) 



If we put Eq. (^T]) into Eq. ( |2B|) and continue in this way for n times, 
we wiU find 



w 



Z = T.[^-{l-q)P{e,-E)]- (32) 

i 

= qt[^-{'^-qme^-E)f^ 

i 

= qi2q-l)t[^-il-q)Pie,-E)]'^ 

i 

= q{2q - l)(3g - 2) ±[1 - (1 - g)/3(e, - E)pl 



q{2q - l){3q - 2)...{nq - n + 1) E[l " (1 " q)P{e^ - E)]^-^ 



with n = 0, 1, 2.... We create in this way a series of equahties which seem 
not to hold. For example, if we take the second equality of Eq. ( |32| ) and 
let g ^ 0, the right-hand side will tend to zero and the left-hand side to 
T,f I = w. The result is w ^ 0. This same result can also be obtained 
for g ^ I if we take the third equahty of Eq. ( |321) and for g ^ | with 
the forth equality and so on. These singular points in q value do not 
conform with the hypothesis that Eq. ( PTp is a basic relation of the 
theory. It seems to us that these equalities are valid only when q ^ 1 
and Z becomes the BGS partition function. 

Now let us suppose a continuous quantity x within < a: < oo replac- 
ing energy in the invariant distribution. Then Z may be given by the 
following integrations 



Z = 



roo 1 

[l-{l-q)/3{x-x)]—dx. (33) 



or 
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z = 



rOO q 

[l-{l-q)l3{x-x)]—^dx. (34) 



In this case, we should put g > 1 for Z to be calculated when x is large. 
The integration of Eq.(^) is always finite. But Eq.(|3^) needs q < 2 
to be finite. If q > 2, the Z of Eq.(|33|) can be calculated while that of 
Eq.(|34|) diverges. This paradox naturally disappears for g ^ 1. 

3. The third problem concerns the calculation of expectation with Eq.(^) 
and the concomitant distribution Eq.(^. Usually, this calculation allows 
to establish a relation between, e.g., internal energy E and micro-state 
energies e, through an E - Z relationship {U = -f^lnZ within BGS). 
This is the crucial step in the statistical interpretation of thermody- 
namics. From Eq.(l^), it is obvious that E calculation is self-referential 
: E = 0[pi{E)] where C[.] is certain mathematical operation. This 
calculus does not make sense for E because Eq.(0) is itself a relative 
distribution with respect to the internal energy E. So E may be arbi- 
trary for a given spectrum e^. To see this, we introduce the distribution 
function Eq.(0) into the expectation Eq.(g) to obtain : 

E = (35) 

Zl-q ^ ^ 

= ^Ee.[l-(l-^Mez-£^)]^ 
^ i 

= - (1 - qme^ - E)]^. - E^[l - (1 - qMe, - E)]^.}. 

Considering Eqs.(^ and (|2T|) , we get 

1 f)y 



Zd(5 

which leads to, instead of the expected U — Z relation 

dZ 



d(3 



= 0. (37) 



So no relation between macroscopic quantities and the correspondent mi- 
croscopic quantities can be found from the expectation definition Eq.(g) 



10 



and arbitrary expectation is possible for any given states. The thermo- 
dynamics connection of this statistical mechanics becomes questionable 
even impossible. 

In addition, Eg. (|37|) gives rise to another problem similar to the second 
one discussed above. Eg . (|T7|) can be easily verified if we take the stan- 
dard Z given by Eq.(||). But If we take the Z given by Eq.(^), we get 
the following relation 

f:(e. - E)\\ - (1 - g)/3(e, - E)p^ = (38) 



or 



2g-l 



If we repeat the same reasoning with the Z of Eq.([20|), we get 



E = ^S^f (40) 



Fi 

^ 3(7-2 



Y.fp 

"' ~ y-w nq-n+l " 
Pi 

which means Y.Y{^i - E) = or E" = E^ez/w if <? = ^ with 
n = 1,2,3.... i.e., we are led to the microcanonical case. This seems 
contradictory to the generality of Eq . (pT]) . 

Above mathematical difficulties or peculiarities may lead to other difficul- 
ties for the theory and seriously harm the reliability of NSM with invariant 
distribution Eq.(|5|). The validity of escort probability and the resulted ex- 
pectation in NSM also become questionable because it gives invariant dis- 
tributions with respect to the expectation when applied to either Tsallis' or 
Renyi's entropy 
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V. AN INVARIANCE RECOVERED WITH A GENERALIZED ENERGY SHIFT 

Though Eq.([Tj) is not invariant with conventional energy shift (e^ + C), 
it can be invariant with a different energy shift we refer to as generahzed 
energy translation. 

We remember the product probability problem mentioned at the begin- 
ning of the previous section. The existence of thermodynamic equilibrium 
in an interacting system described by Tsallis entropy needs not only Eq. 
and the factorization of total probability given by Eq.(^), also an energy 
pseudoadditivity similar to Eq . and consistent with the product law of 
probabihty Eq.(|23) |T1,|7|], that is 

e,j{A + B) = e,{A) + e,(5) - (1 - q)(3e,{A)e,{B). (41) 

and Pij[eij{A + B)] = pi[ei{A)]pj[ej{B)]. Eq. (^) can be written in the con- 
ventional addition form with a generalized g-addition i.e. eij{A -\- 
B) = ei{A) +q ej{B). Let g-exponential distribution Eq.([T]) be noted by 

= expq{-f3e^), we get 

P^{e^ +, C) = ^^^:^expq[-P{e, +, C)] (42) 

= -^^^y^^expq{-(3C)expq{-(3e^) = -^expq{-(5e,) 
= Pi{et). 

where Z{C) = expq{—(3C) and Z{e) = Y.iexpq{—(3ei). So the invariance of 
probability is recovered in a more general mathematical and physical con- 
text for systems with complex interactions or chaotic space time. A complete 
physical comprehension of this "nonlinear" energy shift certainly needs fur- 
ther investigation. 

Here we would like to mention that, in this generalized context, we can 
construct a generalized arithmetic or algebra based on the g-exponential and 
its inverse function g-logarithm (ln^(a;) = ^j—^) with, in addition to +q, 
— and -^q (or / q) corresponding respectively to the following arith- 
metic properties : expq{—qx) = l/expq{x), lnq{x Xq y) = lnq{x) + lnq{y) and 
lnq{l/qy) = — lnq{y). This g-algebra is semi-classical due to its relationships 
containing simultaneously classical and g-operations, but it is useful for NSM 
and consistent with the complex physical circumstance. Considering the fact 
that the classical arithmetic is constructed on the basis of a fragmented world 
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with independent parts, we conjecture that generalizations are necessary for 
the real messy world containing correlated, entangled and overlapped parts. 
We have already seen the development of a «:-algebra based on a /^-statistics 
The above g-algebra is also a possible one for a different physical cir- 



cumstance. Detailed discussion of this mathematical scenario will be given 
in another paper of ours. 

VI. CONCLUSION 

We have discussed some problems related to the energy translation invari- 
ance of probability distribution. It is shown that a probability distribution 
can be invariant only in the thermodynamic limit and only if there is no long 
term interaction or correlation and no relativistic effect is considered. So we 
believe that this invariance is not to be imposed to nonextensive statistics for 
complex systems containing long term interaction. Some theoretical peculiar- 
ities of the invariant distribution reveal that the thermodynamic connection 
of the NSM might be disturbed by this invariance. So the normalized expec- 
tation Eq.(g) becomes problematic for NSM. In addition, the distributions 
given by this expectation is numerically proved different from that predicted 
by the law of large numbers []T7|. All these problems show that it is nec- 
essary to reconsider the role played by the escort probability in NSM. It 
was shown that the invariant property of probability distribution can exist 
with a generalized energy shift which is consistent with the factorization of 
joint probability (product law of probability) prescribed by the existence of 
thermodynamic equilibrium in correlated systems. 
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